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Abstract 

A method for cylinder pressure curve analysis is presented 
which may be used for real time control and diagnosis of 
Diesel engines. By subtracting the pressure samples of a firing 
motor from the non firing case, a difference pressure curve 
can be obtained which contains hidden information on internal 
motor conditions. To use this information, a real time data 
reduction has been implemented which generates symptoms 
like the centre of gravity or the maximum pressure from the 
data samples and may be processed by a classification 
algorithm (L g. an artificial neural network). 

With the presented method it becomes possible to 
separate heat transfer problems from combustion or injection 
failures and to adjust motor characteristics properly. The 
validity of our approach is supported by simulation and experi- 
mental data. 

1. Introduction 

For real time optimization of Diesel engines, for driving at 
maximum torque or for fast failure detection a quality control 
of the Diesel combustion process is very important. 

Thus, in this paper a method is presented which 

• allows to obtain characteristics of the combustion process with- 
! out knowledge of the exact combustion heat curve. Cylinder 
j pressure data is required, but computational costs are drasti- 

• cally reduced. 

thennodyitamics 



edc 

\ 



combustion/ mechanics 













vehicle 
dead) 




— 



cylinder pressure Pi 
Fig. 1: the Diesel engine as a system 



A Diesel engine may be modelled by inputs, outputs and 
internal states. Injection mass % and injection angle 0eb are 
input variables while the dynamic torque T^ is the main 
output signal. The engine speed n is determined by the load 
dynamics and thus not an independent variable. Internally, the 
engine can be separated into 4 blocks: the electronic diesei 
control (edc) unit, the dynamics of combustion, the thermody- 
namic block and, finally, the kinematics (geometry). In the 
notation of system theory, the cylinder pressure signals ft 
represent internal states, see fig. 1. 

When connected to a vehicle, a Diesel engine should 
actually be viewed as a static actuator since its internal 
dynamics are neglec table compared to the dynamics of the load 
(rotational dynamics of the drive shaft, etc.). 



2. Analysis of Cylinder Pressure Signals 

The combustion in a Diesel engine is a very complicated, non- 
linear and stochastic process. A real time analysis of the 
combustion heat curve would be a good way to determine the 
motor condition. As a matter of fact, the centre of gravity of 
the combustion heat curve area depends on the maximum 
pressure P^,, the maximum increase of pressure (dp/dO)^ 
and the maximum gas temperature T„. By sampling the 
cylinder pressure and solving for the heat transfer, it is princi- 
pally possible to extract proper combustion quality symptoms. 
To calculate the transformed energy per rotation angle, one 
has to solve eq.(l) by integration, see [5]: 

dQ dV L 1 fdul + Vf3n] dp 



10 



- m., 

(i) 



where 
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dQJdB combustion heat curve [J/° degree] 
dQ^ldB energy transfer to wall [J/°degree] 



B 

U 



fired fuel mass [kg] 

fired fiiel in residual gas [kg] 

inner specific energy [J/kg] 

stoichiometric air requirement [kg/kg] 

net calorific value [J/kg] 

air for combustion ratio [-] 



However, real time conditions are difficult to match for com- 
putation of the combustion heat curve due to the high compu- 
tational costs. 

Due to the strong relation between cylinder pressure 
signal and the combustion heat curve, the next best signal for 
monitoring is the pressure signal. In the ideal non-firing case, 
the pressure signal of an engine is almost symmetric to the top 
dead centre (TDC). Thus, the centre of gravity in the firing 
case lies close to the axis of symmetry. Its angle as one of the 
possible symptoms is close to TDC and does not change too 
much with variation of input variables m^, or 0^. It would be 
possible to determine this centre of gravity very accurately by 
taking the average of many cycles. This, however, can not be 
done in real-time. 

Therefore, in this work another method is presented 
that uses the difference between the pressure curve in the fir- 
ing and in the non-firing case. Unfortunately, the non-firing 
pressure signal is usually not available in a running motor. 
Thus, one has to either measure mis signal in advance and 
store it in memory or one need to estimate the non-firing case 
from the pressure signal before and after firing. This second 
approach seemed more appealing and was used in mis work. 
The difference pressure curve itself is difficult to interpret in 
real time, though. By mapping it on a set of symptoms, an on- 
line data reduction is accomplished. 



The Centre of Gravity of a Function 

The following graph shows how to compute the centre of 
gravity of the difference pressure curve. 




BDC 



TDC 



BDC 



Fig. 2:Centr f gravity of the pressur curve 



Let BDC and TDC be abbreviations f r the Bottom Dead 
Centre and th Top Dead Centre and let Ap<0) = Pfw^W - 



Pwo flriaJ (0). The coordinates (G^rJ of the centre of gravity of 
the difference pressure curve may be computed by the follow- 
ing equations: 

TDC 

f 0 Ap d$ 

0 = flk P) 



-75c- 



f Ap d$ 

BDC 

where 0 C is the crank shaft angle of the centre of gravity and 

TDC 

I Ap 2 d$ 



1 

2 



BDC 

-JSc- 



C3) 



J Ap dB 



is the pressure at this point. These two values are strong char- 
acteristics for the examination of the combustion. For moni- 
toring the efficiency of the engine the generated torque is 
important, too. 



Difference Pressure Signal Characteristics 

Figure 3 shows the geometric arrangement of an idealistic 

motor with the most important forces and angles. 





F 
















I ^ V 




Fig. 3: Forces in a combustion engine 

The only interesting parameter for the following conclusions is 
the cylinder volume V (see [5D- Taking its derivate leads to: 



dV 

~d¥ 



2 



l 2 /T 



COS0 

- A 2 sm 2 * 



The dynamic torque T^ generated by the gas pressure pW » 
given by the following equation: 



A r 



*dya 



dV 
d& 



(5) 



wh re A is th area of the cross section of the cylinder and ft 
is the environmental pressure. If the torque is integrated ver 
one cycle, a value proportional to the mechanic work of this 
cycle is obtained: 
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AW = 



dd 



(6) 



If th difference pressure Ap is multiplied by dVIdd 9 a symp- 
tom is generated which is proportional to the averaged torque: 



<7) 



The right band side cf mis equation is the work which is pro- 
duced during combustion. If this part is divided by the dis- 
placement of the cylinder V fc , the "mean difference pressure 
Apj" is evaluated: 

AW _ r A _ dV 



Aft 



d$ 



dd 



(8) 



Ap, is an average value over one combustion cycle. The result- 
ing torque is directly related to this value and reaches its 
maximum at the same point. Thus, the following three symp- 
toms for quality supervision of the combustion are available: 



O c The centre of gravity of the combustion 
x c The height of the centre of gravity 
Ap { The mean difference pressure 



In m following sections is shown bow these characteristics 
are related to the combustion process. 



3. Influence of injection parameters 

As mentioned above, a Diesel engine has two manipulated 
inputs: injection angle and injection mass. Between injection 
and the beginning of the combustion, there is a time delay 
which should be known. In mis case, the combustion process 
could be controlled by the injection parameters. For example, 
one c uld drive the motor at maximum torque production. For 
every engine speed and every injection mass there exists an 
injection angle for which the transformation of thermal energy 
to mechanical energy is optimal. For this injection angle die 
motor reaches its maYimnm torque. 




Vc vd VI 

Fig. 4: Onistant-Voiurrte-Cycie: P-V-Chart 



V t : displacement f the cylinder 
V e : compression space 



To derive a theoretical system model, it is important to know 
how the symptoms and "height and crank angle of centre of 
gravity" vary with modification of inputs. To do so, let us 
assume mat the rberrnodynamical process in a Diesel engine 
can be modelled by a constant-volume-cycle, see fig. 4. This 
assumption implies some simplifications to the real case: The 
time for the combustion process is infinite small and during 
the compression and expansion of the gas heat transfer doesn't 
occur. 



P3- 



F2-- 




Fig. 5: Constam-Vohime-Cycle: 
cylinder pressure versus crane-shaft angle 

Figure 5 shows the cylinder pressure for an ideal constant 
volume cycle. During the adiabatic phases the cylinder pres- 
sure is given by the following equations: 



1 -*2: 



3—4: 



(9) 



With the energy equation U 3 
at point 3 is given by: 
R 



Uj + ♦ Hen the pressure 



Py - Pi 



(Qz, * B m ) 



ao) 



with a = RlC r (Q a + and the gross caloric 

value of the fueL If the difference pressure between the expan- 
sion phase and the combustion phase is evaluated, one gets a 
function which depends only on the cylinder volume V and V 2 
(the volume at the begmning of the combustion): 



Ap cv.vj 51 -=77 



V* 



(11) 



Because V, is equal to V 3 , one gets the following equati n: 

After integration over V. the height of the centre of gravity is 
only a function of the volume V 2 . 
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i * 



(13) 



Integration may be performed over the volume or the crank- 
shaft-angle. In both cases, the results are the same. If the in- 
tegration is performed over the volume, one needs to distin- 
guish between two cases: 



O heat transfer before TDC 
O heat transfer after TDC 



dv e rvv.v^+vj 
dv e [vv.v c +vj 



Heat transfer after TDC is the easier case. Solving eq. (13) 
leads to: 



2k 



1 1 



1 - 




2*-l 


v. * y„\ 


- 


\ v ' ' 







(14) 



If heat transfer occurs before the top dead centre, the integral 
in eq. (13) has to be split up into two parts: From V 2 to V c 
and from V c to V, = V e + V h . For this case the following 
result is obtained: 



1- 



-2 



1- 



2k - 1 V 2 



(16) 

Note that for V 2 = V c eq. (14) and (15) are equal. This is the 
trivial case: the combustion happens exactly at the top dead 
centre. 

In fig. 6 the height of the centre of gravity v c is 
plotted versus the cylinder volume at which combustion 
occurs. 















V2 


* 1 * 
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\ V/ 
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* 









Fig. 6: Height of centr of gravity versus cylinder volume 

It is obvious that ir c is a constantly decreasing curve- The 
function is nonlinear, especially in the interesting area before 
the TDC. For the crank angle f the centre f gravity, 0 e , 



similar relationships can be derived. The centre of gravity 
moves to later crank angles if combustion occurs later, i. e. it 
increases with a later injection angle. 

The results of these theoretical derivations indicate 
that x c and © c are not dependend on output variables like 
engine speed n and should only vary with input signals (injec- 
tion angle and mass). 

4. Experimental results 

In this section it will be shown that simulation results from 
former section can indeed predict the behaviour of an existing 
motor. The data has taken from a dynamical engine test stand 
with a Volkswagen Rabbit motor (4 Cylinders, 40 KW, 1.6 1 
stroke volume). 

cylinder pressure [box] 
9 




user 



340" 



420* 



500* 



580" 



Fig. 7 : cylinder pressure curves, 2400 rpm, full load (i 
red data) 



cylinder xueuoie [bar] 




w 340- «w SDor sar 
Fig. 8 : difference pressure curve: P,^ - P^a fm=z 

Fig. 7 shows measured pressure curves for a firing and a 
nonflring engine. The offset between the two curves is caused 
by sensor drifting and is taken care off during real time com- 
putation. The resulting difference pressure curve is shown in 
fig. 8. Note that by subtracting the non-firing pressure actually 
the variati n range of the centre of gravity is amplified. 

In fig. 9, the height of the centre of gravity is plotted 
versus the injection angle and the injection mass. As predicted 
in section 3, this symptom increases with a growing injection 
mass and with a decreasing injection angle. 

In constraint to the theoretical predictions of section 
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4, O c has a slightly different behaviour. Fig. 10 gives the 
measured crank angle of the centre of gravity plotted versus 
the injection parameters. 




Fig. 9: height of the centre of gravity 




Fig. 10: crank angle of the centre of gravity 



Crank angle 9 e is given in degrees relative to the top 
dead centre and is a constantly increasing function in direction 
of the injection angle. Engine speed in fig. 10 and fig. 9 was 
2400 rpm. 

From the experiments, it seems that in contrast to the 
theoretical results, the symptoms have a week dependency on 
engine speed n. 

5. C l a ssific a ti on with neural networks 

So fer, these results are rather encouraging. Therefor, with a 
proper classification algorithm it should be possible to map the 
- varying symptoms to a set of feults. The whole process is a 
pattern recognition algorithm, see fig. 1 1 . 

For nonlinear classification, artificial neural nets 
(ANN) have recently been shown to have desirable properties, 
see [2]. One of the most widely used algorithm is the feedfor- 
ward neural net with baclq>ropaganon training. As a first 
attempt, the three symptoms and engine speed n were used as 
input variables for ANN training. Several configurations of 
this net were investigated regarding their potential for monitor- 
ing and classification of combustion conditions. Configuration 
and real time application issues are under investigation. 



pattern 




map of fiiults 



Fig. 11: Principle of pattern recognition 



6. Conclusions 

The difference pressure signal has been used for real time 
monitoring and for extraction of motor characteristics. Neural 
networks seem to be a feasible tool for failure detection in the 
combustion process, current research deals with optimizing its 
usability for on board diagnosis systems. 
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